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SHARP UPPER BOUNDS ON RESONANCES FOR PERTURBATIONS 

OF HYPERBOLIC SPACE 

DAVID BORTHWICK 



Abstract. For certain compactly supported metric and/or potential perturbations of 
the Laplacian on HP+ 1 , wc establish an upper bound on the resonance counting func- 
tion with an explicit constant that depends only on the dimension, the radius of the 
| ' unperturbed region in H n + 1 , and the volume of the metric perturbation. This constant 

^v^j i is shown to be sharp in the case of scattering by a spherical obstacle. 
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^ . 1. Introduction 

ON . 

For conformally compact manifolds that are hyperbolic near infinity, we now have fairly 
good control over the growth of the resonance counting function. Upper and lower bounds 
have been obtained for various cases in [TJ [5] [T^l H5J HH [Ml HI]- In this paper we 
rN \ develop techniques which can provide a sharp constant for the upper bound, and apply 

■ these specifically to the case where the manifold is a compactly supported perturbation 

of the hyperbolic space H n+1 . The techniques are inspired by Stefanov's recent proof of 
sharp upper bounds on the resonance counting function for perturbations of the Euclidean 
Laplacian [25] . 

Let Ao denote the positive Laplacian on H n+1 . We can write the Green's function 
associated to Aq explicitly: if Rq(s) := (Ao + s(n — a)) , then 



Date: November 12, 2009. 

2000 Mathematics Subject Classification. Primary 58J50,35P25; Secondary 47A40. 
Supported in part by NSF grant DMS-0901937. 



1 



2 



BORTHWICK 




where F is the Gauss hypergeometric function and a := cosh 2 (ic?(2;, z')). From this expres- 
sion we quickly deduce that Rq(s) admits an analytic extension to s € C if n is even, and 
a meromorphic extension with poles at s = —k for k = 0, 1, 2, . . . if n is odd. In the latter 
case the multiplicities of the poles are given by 

, . , s , , (k + 1) ...(k + n- 1) 
1.2) mo(-k = 2fc + ip '- \ '-. 

Let 7?.g denote the resonance set for Er™ +1 (empty for n even), with resonances repeated 
according to multiplicity. The associated resonance counting function is defined by 

N (t) :=#{Ce^ : |C-f| <*}■ 

For n odd, an asymptotic for N (t) is easily deduced by integrating (|1.2[) . For later usage, 
we introduce the constant 

f , 2 -,s. n odd, 

^j(O) ._ ) (n+l)! ' 

1 n even. 

The resonance counting function asymptotics for H™ +1 are then summarized by 
(1-3) N (t) ~ B^t n+1 , 

as t — > oo. 

The main result of this paper concerns the resonance counting function Np(t) for P a 
compactly supported perturbation of Ao. To describe the class of perturbations precisely, 
let 

K := B(0;r ) C H™ +1 

for some ro > 0. We assume that (X, g) is a smooth Riemannian manifold, possible with 
boundary, such that for some compact K C X, we have 

(X-K,g) = (X -K {h9o ). 

In other words, (^",.9) agrees with H n+1 near infinity. Note that X is allowed to have a 
more complicated topology than H™ +1 , as illustrated in Figure [TJ 

Let A g denote the Laplacian on (X, g), and V £ Cfi°(X) with supp(V) C K. We then 
define the perturbed operator 

P ■■= A 9 + V, 
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where some self-adjoint boundary condition is imposed if X has a boundary. Since Rq(s) 
functions as a good parametrix for Rp(s) := (P — s(n — s)) _1 near infinity, it is straight- 
forward to prove meromorphic continuation of Rp(s). We can thus define the resonance 
set Tip, with resonances repeated according to multiplicity, and the associated counting 
function 

N P (t) := #{C G Tip : |C-f|<t}- 
The arguments of Cuevas-Vodev [5] and Borthwick pQ are easily extended to show that 

N P (t) =0(t n+1 ). 

Our goal in this paper is to refine this estimate by producing an explicit constant Bp for 
this bound, which is sharp in the sense that Np(t) ~ Bpt n+1 holds in at least some cases. 

As in Stefanov's work [23], such a result requires a slightly regularized version of the 
counting function. The basis of our estimate is the following relative counting formula: 

(1.4) £ N "®- N °® dt = 2j o a ^dt + ±- J[lo g \r^+ae^\de + 0(lo g a), 

with t(s) the relative scattering determinant for P and a{t) the corresponding relative 
scattering phase. This formula holds for a general class of background manifolds (X ,g ) 
and for much more general perturbations; see Proposition 13.21 

From the relative counting formula (|1.4p . the role that the asymptotic (|1.3[) for No(t) will 
play is clear. The contribution from the relative scattering phase a(t) is similarly easy to 
account for, because it satisfies a Weyl-type asymptotic as t — > oo, 

(1.5) a(t) = \B [ ^t n+1 +0{t n ), 



2 



where 



Bp 



r ] : = 2 r^±3) [ yol(K ' 9) ~ vol(K °' 5o) ] • 



It is for this result that we must require smoothness of g and V. In various asymptotically 
hyperbolic settings, the scattering phase asymptotic was established by Guilope-Zworski 
[14] , Guillarmou [9] , and Borthwick p] . By adapting of the arguments from [I] , we can ex- 
tend the result to the class of perturbations considered here, for a general class of background 
manifolds (Xo,<7o)- 

Once we have the scattering phase asymptotic, the final step in estimating the right- 
hand side of (|1.4[) is to study the integral of log \t(s)\ over a half-circle. It is here that we 
specialize to H n+1 as the background space. With a combination of singular value techniques 
and asymptotic analysis of Legendre functions, we produce a bound 

(1.6) f 2 \og\T(^ + ae w )\d9<B^a n+1 +o(a n+1 ), 
2tt J_* 

with 

(1.7) 4 2) :=^/ f f mxeze ;; o)]+ dxd o, 
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where [•]_)_ denotes the positive part and 

H{a 7 r) :— Re 2a log cosh?' + \/l + a 2 sinh 2 rj — alog(a 2 — 1) 
cosh r — \fl + a 2 sinh 2 



(1.8) 

+ log 



cosh r + yl + a 2 sinh 2 \ 



(2) (2) (2) 

The ro-dependence of B p is approximated by B p « c n e nr ° for ro large, so that Bp is 
roughly proportional to voI(Kq, go). 

The estimate (|1.6|) leads directly to our main result: 

Theorem 1.1. For P — A g + V,a compactly supported perturbation of the Laplacian Aq 
on H n+1 as described above, we have 

(1.9) (n + 1) / ^P-dt<B P a n+1 + o(a™ +1 ), 



i) 



w/iere B P := Bi 0) + B p 1] + B ( p } . 

To highlight the dependence of Bp on P, we note that the constant is dimensional, 
Bp 1 ' 1 depends on ro and on (K,g) only through its volume, and Bp^ depends only on ro. 
None of these components depends on V. 

The factor (n+1) is included in the formula (|1.9|) so that an asymptotic result of the 
same form as (|1.9|) would be equivalent to Np(t) ~ Bpt n+ , with the same constant. Note 
that the only missing ingredient needed to establish such an asymptotic result is a lower 
bound of the same form as (|1.6p . 

To demonstrate the sharpness of Theorem ll.il we consider explicitly the case of scattering 
by a spherical obstacle in H™ +1 , for which X = H" +1 — _B(0; r ), and P = A |x with Dirichlet 
boundary conditions on dX. Figure [2] shows a sample resonance set for a spherical obstacle 
in H 2 . 

Theorem 1.2. If P is the Dirichlet Laplacian on W n+1 — B(0;ro), then 

N P (t) - B P t n+1 . 

Figure [3] shows that resonance counting functions Np(t) for spherical obstacles in H 2 
with several values of ro. These graphs are based on exact computation of the resonances. 
The approximate values of the asymptotic constants for these cases are 

B P ' r 



Already at t — 10 we can see that the behavior of Np(t) is consistent with the predictions 
of Theorem O 

The paper is organized as follows. In ^J2]we develop basic spectral results, such as mero- 
morphic continuation of the resolvent, in a very general "black box" perturbation setting. 
In ^J3]we narrow the context somewhat, in order to establish a nice factorization formula for 
the relative scattering determinant, from which the relative counting formula (|1.4j) follows. 
Another application of the factorization is the Poisson summation formula for resonances, 
which leads to (|1.5|) . The process of estimating the scattering determinant begins in 21 
with a formula that expresses this determinant in terms of the Poisson kernel on HP +1 . In 
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Figure 2. Resonances for the spherical obstacle of radius 7-o = l in H 2 . 
All points off the real axis have multiplicity two; on the real axis the mul- 
tiplicities are more complicated. 



800 



400 




Figure 3. Resonance counting functions for spherical obstacles of radius 
r in H 2 . 



$5] we exploit this relation to prove (|1.6p and complete the proof of Theorem 11.11 A few 
explicit spherically symmetric examples arc considered in which contains the proof of 
Theorem 11.21 Finally, the asymptotic analysis of Legendre functions that is needed for <JS] 
and SJH]is developed in the Appendix. 
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during a stay at the Centre International des Rencontres Mathematiques in Luminy. I am 
grateful to both of these institutions for their support, and also thank the Emory ICIS for 
additional travel support for the Luminy meeting. 

2. Black box perturbations 

In geometric scattering theory, the term "black box" refers to a general class of pertur- 
bations of the Euclidean Laplacian in R n introduced by Sjostrand-Zworski [22]. Although 
in standard usage this terminology is specific to the Euclidean setting, the same abstract 
formulation can be adapted to other settings. In this section, we will discuss black box 
perturbations in an asymptotically hyperbolic context. Our goal is to set up the definition 
of resonances by demonstrating meromorphic continuation of the resolvent, and then to 
prove a global estimate of the counting function. It makes sense to do this in a general 
setting, since only minor changes are required to adapt previously published arguments. 
This section essentially amounts to a review of known results. 

An asymptotically hyperbolic metric on (Xq, go) admits, by definition, a compactification 
Xq with boundary defining function p such that (Xq, p 2 go) is a smooth, compact Riemannian 
manifold with boundary and \dp\ p 2 ga = 1 on dXo. We will assume that (Xo,go) is even 
in the sense introduced by Guillarmou [TU]. This means that the Taylor series of p 2 go 
at p = contains only even powers of p. Under this assumption the resolvent i?o( s ) := 
(A ff0 — s(n — s)) _1 admits a meromorphic to s € C, with poles of finite rank [101 |18j . 

Appropriating the terminology from the Euclidean case, we define a class of perturbations 
of A go as follows. Given a compact K C Xq, we consider the Hilbert space 

H = H (BL 2 (X Q -K ,dg ), 

where Hq is some abstract Hilbert space filling in for L 2 (Kq, dgo). On Tt we consider a 
self-adjoint operator P with domain T> C Ti, satisfying the following assumptions: 

(1) T>\x -k C H 2 (Xq — Ko,dgo)- If u e H 2 (X — K ,dgo) and u vanishes near K , 
then 

(2) For ueV, 

{Pu)\x Q -K a = A go (u\x ~K )- 

(3) As a map H — > H, l/f (-P + i)^ 1 is compact. 

Here the notations 1r : u i— > u\k and 1x -k '■ u l— * u\x -K denote the orthogonal 
projections TL — > Ho and H — * L 2 (Xo — Ko, dgo), respectively. 

We will refer to an operator P defined as above as a black box perturbation of A 9o . Given 
that meromorphic continuation of the resolvent is already known for A 9o , it is relatively easy 
to extend this result to P. 

Theorem 2.1. Let (Xo,go) be an even asymptotically hyperbolic manifold and P a black 
box perturbation of A go . The resolvent Rp(s) := (P — s(n — s)) _1 admits for any N a 
meromorphic continuation to Res > — N + || as an operator p N 7i — > p~ N H., with poles of 
finite rank. 

Proof. The resolvent R g {s) serves as a suitable parametrix for Rp{s) near the boundary. 
Let XOj X, Xi € Co°(A") be cutoff functions equal to 1 on K , such that \ — on the support 
of xi and xa = 1 on the support of %. Let K\ := suppxi, and define 

Pi :=P\ Kl , 
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as an operator on TL\ki with Dirichlct boundary conditions imposed on dKi, so that Pi is 
self-adjoint. We can naturally regard Xi(Pi ~ Z )~ 1 X as an operator on TL. 
Then for z such that z ^ a{P) we set 

(2.1) M(s) = xi (Pi - %) _1 X + (1 - Xo)flo(*)(l - X)- 
Then 

(2.2) (P-«(n-«))M(s)=J-L(s), 
where L(s) = Li(sq) + £2(^7 z o) + ^3( s ) with 



Li(z ) 
L 3 (s) 



= -[A 9)X i](Pi -z ) _1 X, 

= (s(n - s) - z )xi(Pi - zo) -1 ^! 

- [A 5 ,xo]Po(s)(l-x)- 



Our goal is to prove that L(s) is compact and then apply the analytic Fredholm theorem. 
Consider first the error term Li(so), which we can write as 

L x (z Q ) = -[A g , X i]lx- K {Pi - zo) _1 X 

By definition, lx Q -K a {Pi — z o) _1 maps TL to T>\x -k an d we have assumed that the latter 
is contained in H 2 (X — K , dgo). Since [A gg , xo] is first order with smooth coefficients whose 
compact support is contained in X — K , we see that [A 9 ,xo] is compact as a mapping 
H 2 (X — K , dg a ) 1— > L 2 (X — K , dgo). Hence Li(z a ) is compact TL — > TL. 

The black box assumption that 1k (P ~ is compact implies that 1k (Pi ~ is 
compact on TL\k x - And the resolvent identity 

(2.3) (P 1 - z)- 1 = (P 1 - i)- 1 [l+(z- i)(Pi - z)- 1 

then shows that £2(5,^0) is compact on TL. Finally, the error term L^(s) has a smooth 
kernel contained in p°° p' s C°°(X x X). This implies that for N > 0, L 3 (s) is a compact 
operator on p N TL for Res > — N + ^. 

After adding the pieces together, these arguments show that L(s) is compact on p N TL for 
Res > — N + Using the self-adjointness of Pi and the standard resolvent estimate, 

we can insure that ||(Pi — ^o) _1 || is small by choosing Imz large. Similarly, we can make 
||P s (s)|| small by choosing s in the first quadrant sufficiently far from the real axis and 
the line Res = Thus for some s,z we have ||£(s)|| < 1, implying that / — L(s) is 
invertible at this point. The analytic Fredholm theorem then applies to define (/ — L(s))^ 1 
meromorphically on p N TL for Res > — N + ^. The claimed result follows from 



Rp(s) = M(s)(I - L(s)Y 

n 
2 • 



because M(s) maps p N H — > p~ N TL for Res > —N + f . □ 



The fact that Rp(s) admits meromorphic continuation as a bounded operator on TL for 
Re s > f (the N = case) implies, as an immediate corollary, that 

(2.5) <r(P) n (-00, is discrete. 
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Theorem 12.11 allows us to define resonances associated to P as the poles of Rp(s), with 
multiplicities given by 

rnp(C) := rank Res^ Rp(s). 

Then Tip is defined to be the set of resonances of P, repeated according to the multiplicities 
mp. The corresponding counting function is 

N P (t) := #{C G Tip : |C-||<0- 

The remaining goal of this section is to establish an order-of-growth estimate for Np(t). This 
requires first of all that (Xo,go) be hyperbolic near infinity, in the sense that sectional cur- 
vatures all equal —1 outside some compact set. (No resonance bounds are currently known 
in the asymptotically hyperbolic case without this extra condition.) Such asymptotically 
hyperbolic manifolds are even in particular. 

We must also make some extra assumptions of P: 

(i) The operator P must be bounded below, so that the set (|2.5p is actually finite. 

(ii) The singular values of the resolvent of the cutoff operator Pi introduced in the proof 
of Theorem 12.11 satisfy a growth estimate, 



(2.6) ^(OPi-^r^ciimzr^-— , 

for some C independent of z and k. 

The natural way to satisfy the growth estimate (|2.6p is to assume that Ti = L 2 (X, dg) 
for some Riemannian manifold (X, g), possibly with boundary, and that P is an elliptic 
self-adjoint pseudodifferential operator of order 2. Then to establish (|2.6|) we can start by 
using the resolvent estimate (|2.4j) to estimate 

Hk((Pi - z)' 1 ) < I lmz\~i Hk{\P x - z\~i). 

Let Ak ± denote the Dirichlet Laplacian on K%. Since \P\ — has order —1, the operator 
{Ak 1 + 1) 5 |-Pl ™ z\~i is zeroth order and thus bounded on L 2 (Ki,dg). Then (|2.6|) follows 
from 

/i fc (|Pi -z\-i) <n k {{A Kl +1)-') |(A Kl + l)'|Pi-^r 

(The fact that C can be chosen independently of z follows from the resolvent estimate (|2.4p .) 

Theorem 2.2. Let (Xo-,go) be a conformally compact manifold, hyperbolic near infinity, 
and P a black box perturbation of A go that satisfies the extra assumptions (i) and (ii). Then 

N P {t)=0(t n+1 ). 

Proof. This is a fairly minor generalization of the upper bound proved by Cuevas-Vodev [5] 
and Borthwick pQ. This is because for those arguments the interior metric enters only in 
the interior parametrix term, i.e., the first term on the right in (12. 1|) . The difficult part of 
the upper bound analysis involves the terms supported near infinity, and this part of the 
argument applies immediately to P by the assumption that P\x -k = A go . 

To apply the argument from Cuevas-Vodev, we need to check some estimates on the 
interior error terms L 1 (zq) and L 2 {s, z ). For the former, the fact that T>\x -k C H 2 (X — 
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K , dgg) implies that Li(z ) is bounded as a map H — > H {K\— Kq, dg). If A/^^^o denotes 
the Dirichlct Laplacian on K\ — Kg, then we can estimate 

(2.7) v y 

where we can use (|2.3p and (12. 4|) to see that we may take C to be independent of zq. For 
the Li{s, zo) term, we first of all note that (|2.4|) implies 

11^,^)11 < c |B(n ,: a) T 

| lm z \ 

By the assumption (|2.6[) we can immediately estimate 

(2.8) Mfc(i2(s , Zo) )< C M7_il_i^ fc -^T. 

|Imz | 2 

For the argument in [5] one needs to set zq — -fN(n — 7JV) for each N such that \s\ < N, 
so the precise dependence of these estimates on s and z is significant. The estimates (|2.7[) 
and (|2.8[) correspond precisely to the interior estimates eq.'s (2.23-4)]. The proof of [SJ 
Prop. 1.2] then gives a bound 



#{Ce^ P : ICI < r, arg(C-f)e [-7r + £ ,7r- £ ]} < 



n+l 



To fill in the missing sector containing the negative real axis, we apply the argument 
from Borthwick 1 . Here the interior parametrix enters only in the proof of [H Lemma 5.2]. 
The required bound is that for some constant a> n, ||i?p(s)|| = 0(1) for Res > a. Since P 
is self-adjoint and bounded below by assumption, this follows from the standard resolvent 
estimate. The proof of [3 Prop. 5.1] then shows that 

#{C G Tip : ICI < r, arg(C + a - n) G [f + e, ^ - e]} < C £ r n+1 . 

The combination of estimates in the two regions gives the global result. □ 



Remark 2.3. Colin Guillarmou has noted a mistake in the original argument from [12] . which 
propagated through the arguments in [5] and pQ. The faulty claim is that one can choose 
a family of cutoffs {x 1 } such that x % = 1 m some neighborhood of dX and also so that, 
in local coordinates isometric to the unit half-disk in H" +1 , \ l factors as (p(x)ip(y) in the 
coordinates (x, y) e E™ x K + . It is not possible to satisfy these assumptions simultaneously. 

Fortunately, this problem is relatively easy to fix. There are two sets of cutoffs used in 
these proofs. (All three proofs use the same construction.) The inner cutoffs must 
form a partition of unity near the boundary, but are not actually required to factor in local 
coordinates. The essential requirement for the inner cutoffs is that their derivatives satisfy 
quasi-analytic estimates (see [TJ eq. (2.6)] for example), and this is easily obtained without 
reference to a factorization. The local factorization assumption is crucial only for outer 
cutoffs {x\} (with Xi = 1 on the support of We may keep this assumption in place 
because the outer cutoffs do not form a partition of unity. 
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3. Relative scattering theory 



For this section we continue to assume, as in Theorem 12.21 a conformally compact back- 
ground manifold (Xo,go) that is hyperbolic near infinity. The restriction ho — p 2 go\ p= o 
defines a Riemannian metric on dX$, whose conformal class is independent of p. Thus 8Xq 
is commonly referred to as the "conformal infinity" of (Xo,go). A black box perturbation 
P shares the same conformal infinity, since P agrees with A 9o outside a compact set. 

The scattering matrices Sp(s) and So(s), associated to P and A go , respectively, are 
pseudodifferential operators on 8Xq defined as in [TBI H] ■ Away from the diagonal, we can 
realize the kernel of the scattering matrix as a boundary limit of the resolvent: 

(3.1) <S*(s; x, x ) = lim (pp')~ s i?*(s; z, z) for x ^ x. 

where * = P or 0. (This relationship can be extended to the diagonal if one is sufficiently 
careful - see [T5].) This connection allows us to see that Sp(s) and So{s) differ by a smoothing 
operator, as follows. By applying Rp(s) to (|2.2[) from the left, we obtain the identity 

R P (s) = M(s) + R P (s)L(s). 

Then taking boundary limits as in (|3.1[) gives the kernel of Sp(s) on the left, while on the 
right we obtain the kernel of Sq(s) as the limit of M(s), plus a smooth contribution from 
the L(s) term. This implies that the relative scattering matrix Sp(s)So(s) is determinant 
class, and we define the relative scattering determinant 

(3.2) t(s) :=detSp{s)S (sy 1 . 

Let -H*(s) denote the Hadamard product over the resonance set 1Z*: 

(3.3) H m (a) := ]J n+l), 
where 

/ z 2 z p 

E(z,p) := (l-z)cxp zH 1 1 

V z p 

Proposition 3.1. Assume that (Xo,go) is conformally compact and hyperbolic near infinity, 
and P is a black box perturbation of A go satisfying the extra assumptions (i) and (ii) from 
<TJ1 The relative scattering determinant admits a factorization 

(3.4) T(s)=e<M H f-* R H ; {S) 

Hp(s) H {n-s) 

where q(s) is a polynomial of degree at most n+ 1. 

Proof. Since structure near infinity is unchanged from (Xo,go), the arguments of Guillar- 
mou [11] relating the resolvent and scattering pole multiplicities apply to P. Thus the proof 
of [U Prop. 7.2] shows that (|3.4|) holds with with q(s) a polynomial of unknown degree. 
To control the degree, we use the fact that P is bounded from below to obtain 

\\R P (s)\\ = 0(1), for Res > a, 

for some a > n. Then the proof of [U Lemma 5.2] gives that 

\i9 P (s)\ <e c "< s >" + \ 
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for Res < a — n with dist(s, —No) > rj. The same estimate applies to °Dq{s). In the formula 

H (s) Hp(n-s) 

the Hadamard products have order n + 1. Thus the $*(s) estimates imply that \q{s)\ < 
C|s|™ +1+ ' 5 in the half-plane Res < a — n, for any S > 0. Since q(s) is already known to be 
polynomial, the degree of q(s) is at most n + 1. □ 



One nice application of Proposition ^. H is a Jensen-type formula connecting the resonance 
counting functions to a contour integral involving the relative scattering determinant. To 
state this we introduce the relative scattering phase of P, defined as 

(7(0 : =^logr(§ + i£), 

with branches of the log chosen so that cr(£) is continuous starting from cr(0) = 0. By the 
properties of the relative scattering matrix, cr(£) is real and cr(— £) = — cr(£)- 

The following relative counting formula is the asymptotically hyperbolic analog of a 
formula developed by Froese [7] for Schrodinger operators in the Euclidean setting. 

Proposition 3.2. Assume that P is a black box perturbation of (Xo,go) as in Proposi- 
tio ns. 1\ As a — ► oo, 

rN P (t)-N 0{t ) dt = 2 r*M dt+ L / f iogKf +ae -)i^ + o(iog a ). 

Jo t Jo t 27r J-f 

Proof. According to Proposition 13.11 for Re(s) > 5, r(s) has zeros when n — s e IZp or 
s G 7?-o and the latter case occurs only if s(n — s) lies in the discrete spectrum of A So . 
Likewise, poles of t(s) for Re s > 5 occur when either n — s E IZo or s e 7£p, the latter only 
if s(n — s) lies in the discrete spectrum of P. All of these are counted with multiplicity of 
course. 

Let r] denote the contour +t exp(i[— tt/2, n/2])) U +it, ^ — it], as shown in Figured] 
Assuming t is not the absolute value of a resonance in 1Z or IZo, we have 

^-I-(s)ds = N P (t) - N (t) - 2d P (t) + 2d (t), 

where is the counting function for the (finite) set 7?.* H (§, oo) (the resonances coming 

from the discrete spectrum). Evaluating the contour integral yields 

1 / r' 1 f t' 

<p — (s) (is = Im — <p — (s) ds 

2m J n t 2n J v t 



[ o-'(0<% + Im— / 2 -(^ +£ e 2e W e d<9 
M*) + ^ J_[t^log\T(%+te ie )\d6. 



Now if we divide by t and integrate, we obtain the claimed formula with remainder given 

by 

„ f a d P (t) - da(t) , nn 

2 ^ dt = O(loga). 

Jo * 

□ 
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Figure 4. The contour rj. 



Our second important application of Proposition 13. II is to establish the Poisson formula, 
which will lead to Weyl-type asymptotics for the relative scattering phase. Define the 
meromorphic function T*(s) by 

(2s - n) 0-tr[i?*(s) - - s)], 

for s (fc Z/2. The connection between T*(s) and the relative scattering determinant estab- 
lished by Patterson-Perry [3U1 Prop. 5.3 and Lemma 6.7] depends only on the structure of 
model neighborhoods near infinity, and so carries over to our case without alteration. This 
yields the following Birman-Krein type formula: 

Proposition 3.3. For s (fc Z/2 we have the meromorphic identity, 

-d s \ogT(s) = T P (s) - T (s). 

By the functional calculus, T* is essentially the Fourier transform of the continuous 

part of the wave 0-trace (see [T| Lemma 8.1] for the precise statement). By Propositions ^. II 
and 13.31 we can write 



Tp(s) - T (s) = <9 s log 



? ( S ) h p( s ) H o( n - s ) 
H P (n - s) H (s) 



Taking the Fourier transform just as in the proof of [TJ Thm. 1.2] then gives a relative 
Poisson formula: 

Theorem 3.4. Assume that P is a black box perturbation of (X , g ) as in Provosition \3.1\ 

The difference of regularized wave traces satisfies 



0-tr 



cos [t^P-\ 

= \Y ^~ 
2 ^ 



0-tr 



l r)\t\ 



cos It 

«-§)!< 



in the sense of distributions on K — {0}. 
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The desired asymptotics of the scattering phase correspond to the big singularity of 
the wave trace at t = 0. This singularity is very much analogous to that worked out 
by Duistermaat-Guillemin [6] in the compact case. The following result was proven for 
Riemann surfaces, possible with internal boundary, by Guillope-Zworski [141 Lemma 6.2] 
and for higher dimensional asymptotically hyperbolic manifolds (without boundary) by 
Joshi-Sa Baretto [17]. 

Let (Xo,go) be a Riemannian manifold that is conformally compact and hyperbolic out- 
side some compact set Kq C Xq (a more restrictive class than asymptotically hyperbolic). 
Then we consider another Riemannian manifold (X, g), possibly with boundary, with com- 
pact K C X such that {X — K,g) = (Xo — Kq, go). Let A g denote the Laplacian on (X, g). 
We may also include a potential V € C yD (X), supported in K. Given this setup we define 
the operator 

P := A 3 + V, 

acting on L 2 (X, dg) with some self-adjoint boundary condition imposed on the internal 
boundary dX. Clearly P is a black box perturbation of A go , and it satisfies assumptions 
(i) and (ii) of $2] by the remark preceding Theorem [ 



Proposition 3.5. Assume that P = A g + V as described above. IfipE C Q yo (M.) has support 
in a sufficiently small neighborhood ofO and ip = 1 in some smaller neighborhood ofO, then 



(3.5) 
where 



/oo °° 



i-2fe 



Z TT 



do 



r(n±i) 



■0-vol(X, 5 ). 



Proof. By finite speed of propagation we can use cutoffs to split the wave trace into internal 
and external pieces: 



0-tr 



cos (ty/P - n 2 / 4 ) = tr cos (ty/Px- n 2 /A 



X. 



0-tr 



cos 



(tVA go -n 2 /i) (1-x) 



The small time behavior of the first term (which is an actual trace) is given by Ivrii's result 
for compact manifolds with boundary [15] . For the exterior term we can apply [17j . □ 



Using Proposition 13.31 and the Fourier transform relationship between T*(£) and the 
wave 0-trace, we can extract from Proposition 13.51 the asymptotic behavior of the relative 
scattering phase, defined as 

cr(0 :=f logr(i + zO, 

with branches of the log chosen so that cr(£) is continuous starting from <r(0) = 0. By the 
properties of the relative scattering matrix, cr(£) is real and a(— £) = — <r(£). 



Corollary 3.6. As £ -> +oo, 



^(0 = 



(477)-^ 

r(2±2) 



[vol(K,g) - vol(if ,5o)] C +1 + 0(C)- 
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The argument to derive Corollary 13.61 from Theorem 13.41 and Proposition 13.51 requires 
almost no change from that given for n — 1 by Guillope-Zworski [HJ Thm. 1.5], so we omit 
the details. Proposition 13. II (and in particular the bound on the order of r(s)) supplies the 
additional information needed to extend their result to n > 1. The leading coefficient is 
initially given by a difference of 0- volumes, and we use (X — K,g) = (X n — Kq, go) to reduce 
this to a difference of the volumes of K and Kq. 

4. POISSON KERNEL FORMULAS 

Since the asymptotics of a(t) are given by Corollary 13. 6[ application of the formula from 
Proposition 13.21 requires only estimation of \t(s)\ in the half-plane Res > S. To facilitate 
this estimation, we need a more explicit realization of r(s) as a Fredholm determinant. 
This realization will involve the Poisson kernel for the background metric (Xo,ga). For the 
moment we assume only that (Xo,go) is an even asymptotically hyperbolic metric. 

The Poisson kernel can be derived from the kernel of the resolvent Ro{s) by the limit 

Eq(s; z, x') := limp' s Ro(s; z, z'), 

p'-*0 

for z G Xq and x' G dXo- This kernel defines the Poisson operator 

E (s) : L 2 (dX ,dh) - p- N L 2 (X ,dg ), 

for Res > — N + ^, where h is the metric induced on dXo by p 2 go- For / G C co {dX) we 
can solve (A So — s(n — s))u = by setting u = Eo(s)f. Moreover, for Re s > \ with s(n — s) 
not in the discrete spectrum of A go , u has a two-part asymptotic expansion as p — > 0, 

(4.1) (2s - n)E (s)f ~ p n - s f + p s S (s)f, 

where Sq(s) is the scattering matrix. This expansion, for general choice of /, uniquely 
determines the scattering matrix via meromorphic continuation. 

The same construction works for Sp(s). In particular, if we manage to find a family of 
solutions of (P — s(n — s))u = such that 

(2s-n) U ~ p n - s f + p s f, 

for / G C°°(dXo) and s in some suitable region, then Sp(s) can be identified as the map 
/ -> /'• 

Lemma 4.1. Suppose that P is a black box perturbation of (Xo,<7o) with support in Kq. 
Let XI1X2 € C^{X) be cutoff functions such that Kq C {xi = 1} and suppxi C {X2 = 1}- 
The relative scattering matrix can be written as the Fredholm determinant 

r(s) =det(l + Q(s)), 

where 

Q(s) := (2s - n)E Q (s) t [A Q ,X2}Rp(s)[A Q ,xi}Eo(n - s). 

Proof. Since all of the operators in question are meromorphic families, we can restrict s to 
some convenient set like Re s = ^ , s ^ § to avoid poles in the proof. 
Given / G C°°(dXo), consider the ansatz 

(4.2) u=(l- X i)E (s)f + u\ 

as a solution of (P — s(n — s))u = 0. Then P(l — xi) = A (l — xi) implies that 

-[A , X i]E (s)f + (P- s(n - s))u' = 0. 
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After applying Rp(s) on the left, we see that (P — s(n — s))u = may be solved by setting 

u' = Rp(s)[A , X i}E (s)f. 
Using the assumption on supports of \i an d Xi-, we can derive 

(A - s(n - s))(l - xaW = -[A , xaW + (1 - Xa)(P - «(» - *))«' 

= -[Ao,X2]«'- 

This is compactly supported, so that Rq(s) may be applied to give 
(1 - X2)v! = -Ro(s)[A ,X2]u' 

= -Ro(s)[A 0> X2]Rp(s)[A ,xi}E (s)f. 
From this we can deduce the asymptotic behavior of u' as p — > 0, 



II: 



p s Eo(s) t [A 0jX 2}Rp(s)[A ,xi}E (s)f. 



Using the definition (|4.2|) of u and the known asymptotic (|4.1| for Eo{s)f, we thus derive 
the expansion 

(2s - n)u ~ p n -f + p s S (s)f - P s (2s - n)E ( S y[Ao,X2}Rp(s)[A ,Xi}Eo(s)f. 
We can rewrite this as 

(4.3) (2s - n)u ~ p n - s f + (1 + Q(s))S (s)f, 

using the identity 

Eq(s) = -E (n - s)So(s), 
which follows immediately from (|4.1[) . From (|4.3|) we read off that 

S P (s) = (l + Q(a))S («), 

and the determinant follows. □ 

In order to use Lemma [4. II to estimate t(s), we need explicit knowledge of the background 
Poisson operator Eq(s). At this point we specialize to (Xo,go) = H" +1 and work out 
formulas for Eq(s; z,x'). In the usual H" coordinates, z = (x,y) £ K" x R + , we can read 
off immediately from that 

2s _!„-b r(s) 



J5 (s;2;,x') = 2- 2s - 1 7 r"2 



y 2 + | a; — x'\ 2 



r(s-f + i) 

However, our application requires that Eq(s; •, •) be written in geodesic polar coordinates 
and then decomposed into spherical harmonics. The easiest way to do this is to rederive 
Eq(s; ■, •) from scratch. 

In geodesic polar coordinates, H™ +1 = R + x S n and the hyperbolic metric is given by 

go = dr 2 + sinh 2 r dui 2 , 

where dio 2 denotes the standard sphere metric on S n . It is thus natural to adopt the 
boundary defining function 

(4.4) p = 2e- r , 

so that h, the metric induced on dXo by p 2 go, is also the standard sphere metric. 
The Laplacian on H n+1 is 

A = — — ^<9r(sinh" r d r ) + —^A Sn . 
smh r smh r 
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The eigenfunctions of Ag™ are spherical harmonics Y™ with 

AsnY™ = 1(1 + n - 1)17". 
Here I = 0, 1, 2, . . . and to = 0, 1, ... , /i n (Z) with 

4.5 /in l := t—[ „ 

n — 1 \ n — 2 

Proposition 4.2. For EP+ 1 i/ie Poisson kernel in geodesic polar coordinates admits an 
expansion 

oo h„(l) 

(4.6) 2ft(»;r,w,u/) = E E 

(=0 m=l 

wrai/i coefficients given by 

ai (s;r) = 2^-V/ 2 r ( ^ + ^ 1} (sinh r)-^P;_y (cosh r), 
where Pjf(z) is the Legendre function. 

Proof. If we expand the Poisson kernel with respect to the spherical harmonic basis as in 
(|4.6p , then the equation (A 9 — s(n — s))Eq(s) = implies the coefficient equations, 

(4.7) — d^ai — n cothr d r ai + ( — — " — — — s(n — s) ] a; = 0. 

V sinh r J 

After a standard change of variables this becomes the Legendre equation. Since the Poisson 
kernel is smooth in the interior, we select the Legendre solutions that are recessive for r — > 0, 
namely 

a ; (s; r) = Ai (s) (sinh r) 2 P „ +1 2 (coshr), 

for some constants Ai(s). 

The constant Ai(s) may be identified from the asymptotic expansion (|4.ip as r — > oo. 
For the coefficients this expansion implies that 

(2s~n) ai (s;r) ~ p n - s + {So(s)},p s , 

with [So(s)];(s) the matrix elements of the scattering matrix S'o(s), which will be diagonal 
in the spherical harmonic basis. 

Using (|4.4[) and the well-known asymptotics of the Legendre P-function, the leading terms 
in our ansatz as r — > oo are 

_ 1 Y(s — 

(sinh r)-^P s _^+r (coshr) ~ 2 s " ^tt- 1 / 2 -^- ^p"-* 

+ F(Z + n- S ) P ' 

from which we deduce 



□ 
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For future reference, note that we can also read off from this construction the (well-known) 
matrix elements of So(s), 



2s r(f-*) r(z 



(4.8) [*(-)], = * r{s _ 



5. Scattering determinant estimates 



In this section we will combine the formula for r(s) from Lemma 14.11 with the explicit 
Fourier coefficients of the Poisson kernel given in Proposition ^. 21 We can then use estimates 
of the Legendre P-function developed in the Appendix to produce an estimate for the |t(s)| 
term in the counting formula from Proposition 13. 21 

Throughout this section, the background metric is restricted to (Xo,go) — H n+1 . We 
assume that P is a black box perturbation of the hyperbolic Laplacian Ao. As in fjl] the 
support of the perturbation is assumed to lie within 

K Q := {r < r }. 

The main result of this section is the following: 

Theorem 5.1. Let P be a black box perturbation of the hyperbolic Laplacian Aq on H n+1 . 
F, we can estimate 



For a £ § 



as a 



oo, where B 



1 

2^ 

(2) , 



log | 



ae l9 )\d6 <B ( p ] a n+1 + o(a 



n+1 



), 



defined by {1. 7| ). 



Before proceeding with the proof, we note that the combination of Proposition l3.2l Corol- 
lary I3.6[ an d Theorem 15. 1( immediately yields the proof of Theorem 11.11 Note also that 



the restriction to P = A g + V (from the more general black box class) is needed only for 



Corollary 13. 61 

The proof of Theorem 1 5 . 1 1 will be broken into several stages, starting with: 

Lemma 5.2. Assuming P and r$ are defined as above, fix some small e > and rj > 
and define rj := ro For Res > 1| with dist(s(n — s),a(P)) > e, the relative scattering 

determinant can be estimated by 



(5.1) 

where 

(5.2) A,(«):=|2«- 



log|r(s)| <J2 M01og(l + C7A Z ( S ; 



|a;(n — s; r)\ 2 (sinhr)" dr 



\ai(s; r)\ 2 (sinhr)™ dr 



with a;(s;r) the coefficients from Proposition ^. 2\ and C depends only on e, r\, and rg. 

Proof. Let xi an d X2 be smooth cutoffs as in Lemma l4.1| such that Xj = 1 f° r r ^ r j an d 

Xj 



for r > rj+i. Then we can rewrite the Q(s) from Lemma 14. II as 

Q(s) = (2s - n)Eo{s) t l lr2n] [Ao,X2}Rp(s)[A ,Xi}Mr 1 ,r 2 ]Eo(n - s), 
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where l[ ri , r <+1 ] denotes the characteristic function X[ri,r i+1 ]( r )i acting as a multiplication 
operator. By Lemma B~T1 and the cyclicity of the trace we have 



log|r(s)| = det(l + (2s - u)[A ,X2]i?p(s)[Ao,xi]l [ri , r2 ]^o(n - s)E (s) i*l [raj1 

For Res > under the assumption dist(s(l — s),cr(P)) > e, we can apply the spectral 
theorem and standard elliptic estimates to obtain 



[A ,x2]Rp(s)[A Q ,xi] 



<c, 



where C depends on e, 77, and ro. Under these restrictions 
(5.3) 

where 



log|r(s)| <5>g(l + ^(F(s))j, 

3=1 



F(s) := (2s - n) l [rur2] E Q (n - s).E (s) l[ r2 ,r 3 ] 



Using Proposition ^. 21 the eigenfunctions of F*F(s) can then be written down explicitly. 
If we define 

ui, m (r,oj) := l [r2!l . 3] a/(s;r)r, m (w), 

then 

F*F(s)ui tm = A;(s) 2 u iim , 
where Xi(s) is given by (|5.2p . To see that {A;(s)}, counted with multiplicities, contains 
all of the nonzero eigenvalues of F*F(s), suppose that w € L 2 (H" +1 ) and (ui, m ,w) = 
for all l,m. Then by (|4.6[) we have £Jo(s)*l[r 2 ,r 3 ] w = 0, which implies that F*F(s)w = 0. 
Hence, after possible rearrangement, the sequences {Aj(s)} and {/ij(F(s))} correspond. The 
claimed estimate follows from (15.31). □ 



Lemma 15.21 reduces our problem to the estimation of the A;(s)'s, which we take up next. 

Lemma 5.3. Assume that Res > § and \s — %\ eN, and set k = I + and ka — s — § . 
Assuming that r% e (ro,ro + 1) in the definition \5.2\l o/A;(s), we have the bound 

logAi(s) < kH(a,r 3 ) + Clogfc, 

where H(a,r) was defined in hi. 8^ . with a constant C that depends only on ro- 

Proof. From Proposition 14. 21 we obtain the explicit formula 



Aj(a) = sin7r(s - f ) T(l + s)T(l + n - s) 



P~ 



(cosh r) 



sinh r dr 



t n-l 2 

P s , 1+ i 2 (coshr) 



''2 



sinh r dr 



where we have exploited the symmetry P~ k {z) = PZx- v { z )- 

By conjugation, if necessary, we can assume that arga <E [0, f ]. Applying the estimate 
from Corollary IA.2I then yields 

sin(Trfca) T(fc(l + a) + \) T{k{\ - a) + \) 



(5.4) 



r(fc + i) 2 

e 2fcRo(0(a,r)-p(a)) sinhr dr 



s 2feRe(0(a,r)-p( a )) ^ dj . 
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By (|A.6[) . Re<fr(a,r) is increasing as a function of r. Hence 



s 2fcRe(*(a,r)-p(a)) sinhr dr 



< e feRcWa ' r3) - p(Q)) coshr 3 , j = 1,2. 



The first factor on the right side of (|5.4p can be estimated directly via Stirling's formula 
for a £ [1, oo), 

sinvrfca r(fc(l + a) + ±) T{k(l - a) + |) 



(5.5) 



log 



r(fc + i) 2 



7rfc| Imal + fcRe 



(a + 1) bg(a + 1) + (1 - a) log(l - a) + 0(log fc), 



as A: — > oo, uniformly for arg(a— 1) > 5. We can extend the same estimate to arg(a — 1) < 8, 
using 



sin7rfca — a) + 5) 



— TTtsainka 
T(k(a-l) + ^) , 



and our assumption that |fca| S N, which implies 

I tan7rfca| < 1. 

After we note that 

H(a, r) = Re 2(j>{a, r) - 2p(a) + (a + 1) log(a + 1) - (a - 1) log(a - 1) 

we obtain from (|5.4[) and (|5.5[) the estimate 

logAz(s) < kH(a,r 3 ) + O(logfc) + 2 log cosh r 3 . 



□ 



Now we can combine Lemmas 15.21 and 15.31 to estimate t(s). The strategy here is similar 
to Stefanov's in [23l Thm. 5a]. 

Proposition 5.4. For a — 5 6 N and l^l < | we ftai>e 

log |r(f + ae 4e )| < 6(0, r )a n+1 + o(a" +1 ), 
uniformly for \9\ < ^ — ea~ 2 , urai/i 



(5.6) 



r(n) Jo 



[g(xe'V )]_ 

~.n+2 



where [-] + denotes the positive part. 
Proof. Since 



s(n — s) 



n 

T 



a 2 e 2ie 



the assumption that \9\ < 5 — ea 2 implies that (s(n — s) remains a distance 0{e) from <r(P) 



for a sufficiently large. The hypothesis of Lemma [572] is thus satisfied, yielding the estimate 
(|5.ip with a C that depends only on e, 77, and ro. To apply Lemma 15.31 to estimate the 
right-hand side of (|5.1|) . we need to distinguish the terms according to the sign of H(a, r 3 ). 
For large a the sum is dominated by terms with H(a, r 3 ) > 0, which occurs for a outside a 
certain neighborhood of the origin, as shown in Figure GO 

Let x = A(9) be the implicit solution of the equation H (xe 10 , r 3 ) = 0, so that H (xe lS , r 3 ) > 
precisely when x > A(9). Given some <5 > 0, we will subdivide the sum (|5.ip by breaking 
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Figure 5. The positive region for H(a,r), shown for r = 1. 



at values where \a\ = A(0) and (1 — 5)A{6), leaving us with three parts. The dominant part 
of the sum will be 

£+:= K(l)log(l + CXi(s) 

I: \a\>A(8) 

(Recall that a = (s — §)/A; where k = I + (n — l)/2.) For a in this range, assuming 
I arga| < ~, we apply Lemma 15.31 to obtain 

(5.7) log(l + C*Ai(s)) < kH(a,r 3 ) + Clog k 
Using this estimate together with the asymptotic 

om-1 

T(n) 

we have 

(5.8) s+ < y: (w +CFl ")( fci/ (T-' r3 ) +ciogfc 

k<a/A(9) V V ' / V 

(The sum could be restricted to k > ^y^, but this would not improve the bound.) We can 
estimate H(a, 7-3) = 0(\a\) with a constant that depends only on r%. Thus 

k<a/A(9) 

With this estimate, the sums over lower order terms in (|5.8p are easily controlled, and we 
obtain 

(5.9) E+ <_JL £ k n H (^^ +Ca n logaj 

V ' k<a/A{9) 

where C depends only on e, rj, and tq. 
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Because H(xe > , r) is an increasing function of x, the right-hand side of (|5.9[) is easily 
estimated by the corresponding integral, 

2 /"afar 

£+ < 



r(n) Jo 

Making the substitution x = a/k gives 



r A{6) / np™ \ 

J k n H[^ — ,r 3 jdk + Ca n \og> 



2a"+! f° [g(se« r 3 )]+ J 

E+ -Wio ^ + g ° 

= &(0,r 3 )a n+1 + Ca™ log a, 



with C depending only on e, 77, and ro. 
The middle term in (15.11) will be 



E := ^ h n (01og(l + CAKs)). 

/: (l-8)A(8)<\a\<A(8) 

The number of terms in this sum is 0{a8), and we can control them using (|5.7p . noting also 
that H(a,r 3 ) = 0(5) for |a| in the given range. Using an integral estimate as we did for 
we thus obtain 

S < C5a n+1 +Ca n logo, 

where C depends only on e, tj, and ro- 
The final portion of the sum is 



Y, Mi)iog(i + CA,00 



£_ := 

(: \a\<(l-S)A(9) 

We use the fact that H{a, r 3 ) < —CS in this range to estimate 

log(l + CA;(s)) < CX^s) < Ce- Ck . 

This implies 

£- < C s e' ca , 

for some c > 0, where C$ depends on 5 as well as e and the r,-'s. 
Adding the three parts S+, £ , £_ of (|5.ip together now yields 

(5.10) log|r(| +ae l6 )\ < b(6, r 3 )a n+1 + C(e, 77, r ) [<5a" +1 + a n log a] + C(e, 77, r , <5)e~ ca , 

where we have made the dependence of the constants explicit. Since H(a 1 r) is a strictly 
increasing function of r, we can absorb the Sa n+1 term into the first term by assuming that 
5 is small relative to 77 and replacing r 3 with = ro + 4r). With this change, we obtain 
from (|5.10p the estimate 

log \t(- + ae i9 )\ 

(5.11) S ' V2 „ +1 < b(6, r + 47,) + C(e, rj, r , ^a" 1 log a. 

The constant C(e, 77, ro, 6) may well blow up as 77 — > 0. The best we can do here is to observe 
that (|5.1ip implies 



lim sup 



log|r(f + ae^)| 

^T+i 6 ^' r °) 



<6(0,ro+ 477)-6(^,r o ), 



for any 77 > 0. Since b(9, r) is uniformly continuous on [— § , -|] X [ro, ro + 1], we can now let 
77 — > to obtain the claimed o(a" +1 ) estimate. □ 
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Proof of Theorem \5.1\ For any e > 0, we can integrate the result from Proposition [52] over 
\9\ < ^ — ea~ 2 , which gives, 

/ log |r(f + ae t6 )\ dS < B [ p ] a n+1 + o(a n+1 ). 



The factorization given by Proposition 13.11 together with the minimum modulus theorem 
(see e.g. [24] Thm. 8.71]), implies that for any <5 > 0, there exists a sequence fj — * oo such 
that 



(5.12) |r(f + ne ie )\ < C s e X p(r?+ 1+5 ), 

uniformly in 9. In sectors of the form \6\ £ [§ — /3, §], where + ae l9 ) is analytic, we 
can apply a Phragmen-Lindelof argument, using (|5.12[> . log \t(s)\ — for Res = §, and the 
estimate from Proposition 15.41 for |0| = 5 — f3, to conclude that 

\ T { r j + ae l9 )\ <Ca n+ \ 
uniformly for |0| G [f - j3, § ]. Thus, 

/ log |r(f + ae t9 )\d8< Cea n ~\ 

-ea- 2 <|6»|<f 

□ 



6. Examples 

Suppose that X = M n+1 and we consider a black box perturbation P = A g + V, where 
both the metric g and potential V are spherically symmetric. The symmetry assumption 
guarantees that the perturbed Poisson kernel is "diagonalized" by spherical harmonics, in 
the sense that 

oo MO 

1=0 m=l 

The coefficients a/(s;r) will satisfy (|4.5[) for r > tq and are thus expressible in terms of 
Legendre functions. Following the convention of Olver, we use the Legendre Q-function in 
the form 



(6.1) 



Q2(*) 



r(v + fi + i) 



where Q£(z) is the standard definition. This makes Q^(z) an entire function of either /i or 
v, which is much more convenient for identifying resonances. We can formulate the general 
solution of (|4.5[) for r > r as 



(6.2) 

where 
(6.3) 



ai(s; r) = (sinhr) 



4(s)Q*(ooshr) + £,(s)Q* ^(coshr) 



/ + 



n-1 



n + 1 
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In particular examples, Ai(s) and Bi(s) will be determined by matching a; and its first 
derivative to the corresponding solutions for r < tq . The scattering matrix elements can be 
read off from the asymptotics of these solutions as r — > oo, using 

o/(a;r)~e ! .(p n - + p > [5p(«)] J ) > 

in the same way that we found [So(s)]z in (|4.8[) , Indeed, from the well-known asymptotic 
[H eq. (12.09)] 

1 —u—l 

(6.4) Qfo) = *' + 3 (|) (l + 0(z- 2 )), as*-oo, 
we can see from (|6.2p that the scattering matrix elements are given by 

(6.5) [Sp(s)], = 



t _ 2s r(f -a) Ms) 



Consider the case where P is the Laplacian for a spherical obstacle of radius r$ in H ra+1 . 
Imposing the Dirichlet boundary condition at r = ro gives coefficients 

A s = Ql^coshro), B s = -Q^(coshr ). 

In this case, from (|6.5p we see that 



(6.6) [S P {a)]i 



, n -2s r(§-a) QVi(coshro) 
J r(s-f) Q^(coshro) ■ 
With this observation we can give the: 

Proof of Theorem \1.2l Our goal is to show that 

(6.7) ±- / 2 log|r(f + ae ie )\d9~B n (r )a n+1 . 

In conjunction with Proposition 13.21 and Corollary 13. 61 (|6.7I) would imply that 

(„ + l) f a ^l dt ^B P a n+1 , 
Jo t 

and this is equivalent to the stated asymptotic for Np(t). 

Using (|6.6p with (|4.8p gives the relative scattering matrix elements, 

[ ^ (S)50(S) Q^(coshro) ' 

with k and ^ defined as in (|6.3p . With the connection formula [TH1 eq. (12.12)], 

Q*„-i(z) Qg(g) , v p _„, v 

— 71 7 T, T = COS (-KV) P., A (Z , 

r(fc + ^ + i) r(fe-i/) v > v \ h 

we can rewrite the coefficient in the form 

(6.8) [S P (s)S (s)- 1 ] l = 1 - cos(7ri/) T(k v /^J™^ ^ ■ 

Q£(cosh r ) 

Now consider 



log|r(a)| = 5^(0 l°g [M^oW -1 ] 



1=0 
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Defining a by ka — s — tt , we can use (|6.8p to write this as 

oo 

(6.9) log |r(s)|=£ /»»(!) log |1 -m(a)l 

where 



z=o 



rjk{oi) :— sin(7rfca) T(k(l — a) + g) 



^ri +fca (coshr ) 
QN +fca (coshr ) 



Assuming that | arga| < -| — e, Corollarv lA.3l gives the estimate 



log | Vk (a) | 



sin(Trfca) T(fc(l - a) + ±)r(fca + 1) 



= fc Rc[20(a,r o )-p(a)— <?(a)] 



r(fc + i) 

with constants depending only on s. Applying Stirling's formula and avoiding the poles by 
assuming | s — ^ | £ N as in the proof of Lemma 15.31 we have 

sin(vrfca) T(k(l - a) + ±)T{ka + 1) 



r(* + i) 



fcRe 



a log a — (a — l)log(a — 1) +0(loga) 



Since 



2(j)(a, r ) — p(a) — q(a) + a log a — (a — 1) log(a - 1) = H (a, r ), 



log|77 fe (a)| x kH(a,r Q ) + O(loga). 



Re 

the full estimate is 
(6.10) 

As in the proof of Proposition 15.41 we divide the sum (|6.9[) into three pieces £+, Eo, and 
with breaks at |a| = (1 ± 5)A(Q) for some S > 0. The dominant piece is 

£+:= h n (l)log\l-r]k(a)\ 

l: \a\>{l+5)A(e) 

Using the lower bound from (|6.10p . but otherwise arguing as in the proof of Proposition l5.41 
we have 

£+ > b{d,r Q )a n+1 +0(a n log a) 
The estimates on Eo and E_ are identical to those in Proposition [5]4] 

£ < cSa n+1 +0(a n logo), 

and 

£_ = 0{e~ ca ). 

Hence we conclude that 

log|r(f + ae ie )\ > (b(0, r ) - cS)a n+1 + 0(a n logo), 

for a £ N and \8\ < \ — e, with constants that depend only on ro and e. 

Integrating, over 9, and using Proposition l5.4l to control the errors from \9\ £ [f — £, § ], 
we obtain the estimate 



1 



— / log |r(f + ae w )\ d6 > (B^ - e)a n+1 - C nue a n logo, 

valid for any e > 0. (This e combines the terms proportional to e and S from above.) In 
combination with Theorem 15. 1[ this proves ()6.7|) . □ 
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With some care, the explicit scattering matrix provided by (|6.5p can be used to compute 
resonances. Scattering poles and zeros are defined a renormalized scattering matrix Sp(s), 
in which the infinite rank poles and zeros coming from the gamma functions are removed, 

The scattering multiplicity is then defined by 

MO ■= -tT[Res c S'p^Spis)- 1 ]. 
For H™ +1 , the connection between scattering multiplicities and resonances is given by [HJ 

canu 

n odd 

K/2-1 (C) - K/2+1 (C)) h n +i(l) n even 

For Re£ < ^, the term mp{n — Q plays a role only if P has discrete spectrum. In the 
examples that we will consider explicitly, n = 1 and the discrete spectrum is empty. For 
these cases, the resonances are precisely the poles of the Sp(s). 

Consider first the spherical obstacle of radius r in H 2 , for which the scattering matrix 
is given by (|6.6|) . From this expression we can read off the resonance set 

K P = (J{s:Qti(coshr ) = 0}. 

fcez 

Figures [2] and [3] were thus obtained through numerical computation of zeroes of the Legendre 
Q-function. 

As a second example, we consider scattering in H 2 by a radial step potential of the form 



v p{Q = m p(C) - m P (n - C) + 




V(r) = 



c r<r , 
r > tq. 



In this case, with P = Aq + V, the coefficient solutions for r < rg are Legendre P functions 
PjA(ro), with 



oj(s) :=-I + ^( s -I)2 + c . 
The corresponding resonance set is 

(6-11) Kp=\j{s:W[Q k ;_ 1 (z),p- k (z)] = o}. 

I K ' 2-coshrn ' 

where W is the Wronskian. Resonance counting functions for c = 1 and c = 5, with ro = 1, 
are shown in Figured] 

We should note that Theorem 1 1 . 1 1 does not apply to the step potential, because the lack of 
smoothness means that we cannot derive scattering phase asymptotics through Corollarv l3.6l 
In view of the scattering phase asymptotics proved by Christiansen [4] in the black box 
Euclidean case, one might hope that the smoothness requirement in our case could be 
loosened. However, the technique of Robert used in [4] does not seem to be applicable 
to the conformally compact hyperbolic case. In any case, it is interesting to compare the 
putative upper bound suggested by Theorem 11.11 to the empirical results based on (|6.11[) . 
For both of the cases shown in Figure [6l the constant from the theorem would be Bp « 
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Figure 6. Resonance counting functions for radial step potentials in H 2 . 



3.15. The numerical results thus suggest that Np(t) satisfies an asymptotic with a constant 
significantly smaller than upper bound that Theorem 11.11 would predict. 

Our final example is a "transparent" spherical obstacle. Let P = A g where 

)K 2 g r<r , 



Then 



where 



[.9o r > r Q . 
Tip = \j{s : W[Q k s „i(ro),P- { k s) (rv)] = o}, 



cu(s) : =-I + v / K 2 s(s _ 1)+ i 

Figure [7] shows resonance counting functions for n — i and n = Once again, Thco- 
rcm ll.ll does not apply because of the lack of smoothness. However, in this case the predicted 
constants, 



B P = 



2.75 k=|, 
3.70 k = V2, 



at least roughly match the observed behavior, so that one might believe that the theorem 
would give a sharp result if extended to this case. 



Appendix A. Legendre function estimates 

In this section we will estimate the growth of the Legendre functions P^(coshr) and 
Q^(coshr) as k, \v\ — > oo simultaneously. We wish to extract the leading asymptotic behav- 
ior, with error bounds uniform in a := (y + \ )/k for Re a > 0. The construction of these 
estimates leans heavily on techniques from Olver [19| §11]. 
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Figure 7. Resonance counting functions for transparent spherical obsta- 
cles in H 2 with ro = 1. 



Throughout this discussion we identify z — coshr and switch freely between the two 
variables. Let 



w{z) = (sinhr) 



P_L fc (coshr), 



Q"| +fcQ (coshr) 



Then the Legendre equation reduces to 

(A.l) d 2 z w=(k 2 f + g)w, 

with 

1 + a 2 sinh 2 r , „ sinh 2 r 



(A.2) f { r)^—— 4 , g{r):=- . 

smn r 4smn r 

If Re a = then the equation (|A.1[) has turning points (points where / vanishes to first 

order) when a = ±z/ sinhr. By conjugation, it suffices to assume Imct > and so we focus 

on the upper turning point. To obtain uniform estimates near this point, we introduce the 

complex variable £ defined by integrating 

(A.3) ^dC = y/fdz, 

starting from ( = on the left and from zq = — I /a 2 (the turning point) on the right. 
Throughout this section we assume principal branches for the logs and square roots, under 
the restriction that arga G [0, 7r/2]. 
Integrating both sides of (1A.3|) yields 

(A.4) §C« = 0, 

where 

,. . f r y/l + a 2 sinh 2 1 , + 

(p{a,r) := / at 

J cosh _ 1 



(A.5) 



sinh£ 



/ a coshr + \J 1 + a 2 sinh 2 r \ 1 , 

alog 7^=i + 2 log 



cosh r — Vl + ot 2 sinh 2 r 
cosh r + y/ 1 + a 2 sinh 2 r 
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Figure 8. Trajectories of <fi(a, •) with |a| = 2 and 9 = arga. The turning 
point occurs at the origin. 




Figure 9. Trajectories of ((a, •) with \a\ = 2 and 9 = arga. 

The expression (|A.5[) is well-defined by principal branches for arga € (0, 7r/2], and we extend 
the definition to the positive real axis by continuity. (At the apparent singularity at a = f , 
this extension yields </>(f,r) = logsinhr.) 

The region of interest, namely arga € [0,7r/2] and r > 0, corresponds to the sector 
arg^> e [— 7r, as illustrated in Figured] Figure [9] show the corresponding picture for £, 
and illustrates in particular how passing from <p to C resolves the singularity at the turning 
point. For future reference, we note that <j) satisfies the equation 

(A. 6) d r (j> = \fj sinhr, 

implying in particular that Re d r 4> > 0. The fact that Re <p is an increasing function of r 
will be important later, and is not so evident from (|A.5[) . 
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The asymptotics of <fi(a, •) will also play a crucial role. As r — + 0, we have 
(A.7) <f>(a,r) = log (~) + p(a) + 0(r 2 ) 1 

where 

(A.8) p(a):=|log 
And as r — > oo, we have 



a + 1 
a- 1 



•log(l-a 2 ), 



(A.9) 
where 

(A.10) 



(a, r) = ar + q(a) + 0(r 2 ), 



q(a) := a log 



log 



l-a 
1 + a 



Proposition A.l. Assuming that k > 0, arga €E [0, 5] and r € [0, oo), we have 



i 27r 2 fc e (T 4 e 8 

<*■"> ^'- h ^ r (t + 1 ) [1 + a : 5i „ tfr]t 

and 

(A.12) Q*, (coshr) 



-kp(a) 



Ai(/c3e 2 "C) + hi(k,a,r) 



2n 



r(fc« + i) [i + a 2 sinh 2 r ]i 



kq{a) 



Ai(/c^C) + h (k,a,r) 



(A.13) 



where £ is defined by $A.4\) and KA.5\) , p(a) and q(a) are defined in iA.fy) and HA. 10\) , 
respectively. The error terms satisfy 

\k*Cihi(k,a,r)\ < Ce^^Ar 1 (l + |a| - *V 

\ki(lho(k,a,r)\ < Ce^^hT 1 (l + |a|-l) . 
ure'i/i C independent of both a and r. 

Proof. If we set W = (f / (.Y^w, then the equation (|A.1[) transforms to: 
(A.14) 9 2 iy = (fc 2 C + 

a perturbed version of the Airy equation, with the extra term given by 

c 



(A.15) 



Following Olver P3B Thm. 11.9.1], we consider solutions of the form 
(A.16) W a = A\(kie^Q + h a (k,a,r), 

for a = — 1, 0, 1, where the error terms satisfy 



(A.17) 



djh a - k 2 (h a = ip h a + Ai(fc3 e ^c) 



Let us focus first on the Legendre P-function. As \w\ — ► oo, the Airy function Ai(w) is 
exponentially decreasing for |argw| < ? and exponentially increasing for |argu;| € (-|,7r]- 
Since -P~* +fca (cosh r) is recessive at zero, and r — > corresponds to £ — > e^^oo, we choose 
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the solution W± from (|A. 16|) . The assumption that the solution is recessive as r — > implies 
boundary conditions. 

(A.18) Jn|r=o=$-Jn|r=o=0, 

which we must impose on (|A.17[) . To identify the Legendre P function with a multiple of 
W\ , we compare the well-known asymptotic 

p-k (coshr)^^i — (l\ k (i + 0{r 2 )), 

to the behavior of the ansatz 

(sinhr)- 1 ^//)^! = \ (Ai(fc§e 3 ^C) + fti(fc,«,r))- 

Away from the negative real axis, the Airy function has the asymptotic behavior [19, 
eq. (4.4.03)] 

(A.19) Ai(u>) = -t-rw-i cxp(-fu>i) [l + 0(\w\-§)] , 

27T2 

uniformly for | argw| < 7r — 5, with a constant that depends only on 6 > 0. To cover the 
negative real axis we have also [TH1 eq. (4.4.05)], 

(A.20) Ai(«;) = X(-w)-i cos(f (-«,)! - f ) [1 + 0(M-*)] , 

uniformly for | argw| S + 5, 7r]. (These estimates agree where they overlap.) 

As r — > 0, we have — > +oo, which is in the range covered by (|A.19|) . Along with 
the asymptotic behavior of £ deduced from (|A.7[) . this yields 

(iAi(kh^() ~ fc-s^U fcp(Q) (-) as r -> 0. 

27T2 \2/ 

Thus for the Legendre P-function we find 

^rJ +te (coshr) = 27r^i e y r - (sinhr)- 1 ^//) Vi, 
which proves (|A.11|) . 

To complete the analysis of the P case, it remains to control the size of the error term 
hi(k, a,r). The error bounds may be derived as in the proof of [19[ Thm. 11.9.1], starting 
from the differential equation (|A.17p satisfied by h\. Using the boundary condition (|A. 18|1 . 
we can apply variation of parameters to transform this to an integral equation, 

h x {k, a, r) = -^£1 f Kl (r, rW) [hx(k, a, r 1 ) + Ai(fc§ e*P((r'))] ^ smW dr > , 

K3 Jo C(' ) 2 

where 

Kx(r, r') := Ai(kh^C(r'))Ai(ki((r)) - Ai(ki<;(r'))Ai(kh^((r)). 

Then, using the method of successive approximations as in pj|l Thm. 6.10.2], together with 
the bounds on the Airy function and its derivatives developed in pjU §11.8], we obtain the 
bound. 



(A.21) fc*C*fti < Ce feRc<#, (e cfe_1 * lW - l) 
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where 
(A.22) 



Jo 



sinh r . 



Using (|A.2I) and (|A.15|) . direct computation shows that 



(A.23) 







- sinh 2 



4a 2 cosh 



4(1 + a 2 sinh r) 2 



5 (1 + a 2 sinh 2 r) 2 



1G 



£2 sinh r 



Because we require estimates that are uniform in both a and r, the analysis of (|A.23[) is 
somewhat complicated. For some small c > 0, we will break the estimation into 3 different 
zones as described below. We use the notation A x B to mean that the ratio A/B is 
bounded above and below by positive constants that do not depend on a or r. 

Zone 1: Assume that |1 + a 2 sinh 2 r| > c and |a| > 1. The first term in the formula (|A.5|) 
for <j> dominates for large r and the second term for small r. We can thus derive the bounds, 



\4\ 



log|a|r for |a| sinhr < | 
i\r for \a\ sinhr > |, 



Using this to estimate C = (§0)^ m (|A.23|) gives 



sinhr < 



Ci|a| 2 (— log \a\r) 3 r + Ci{— log |a|r) ^r 1 for |a| sinhr < |, 
Ci|a|~^r^e~ 2r + C2|a|~ir^^ 



for |a| sinhr > i. 



It is then relatively straightforward to control the contribution of these terms to (|A.22[) . For 
lal > 1, we obtain 



(A.24) / W'C" 

J |1+Q 2 sinh 2 r|>c 

for some C that depends on c but not on a. 
Zone 2: Assume that |1 + a 2 sinh 2 r| > c and 



sinh r dr < C, 




a\ < 1. In this case we claim that 

log(l — e~ r )| for |a| sinhr < i 
a\[r + log2|a|) for |a| sinhr > |, 

Using these in conjunction with (|A.23j) then gives 

C\\ log(l — e~ r )\i sinhr + Ci\ log(l — e -r )| - ^ (sinhr) -1 for \a\ sinhr < 
Ci|a|-§(r + log2|a|)^e- 2r + C 2 |a| - i(r + log2|a|)~" r 
For I a I < 1, we obtain 



sinhr < 



for lal sinhr > 



(A.25) 



/ 

Jll+a 2 



sinhr dr < C\a\ 3 . 



1>f*C 

' 1 1+a 2 sinh 2 r| >c 

for some C that depends on c but not on a. 

Zone 3: Assume that |1 + a 2 sinh 2 r| < c. This puts us near the turning point. It is 
convenient to use the z — coshr variable here. The turning point occurs at the point 



z a := y/l-a 2 , 

which lies near the path of integration for (|A.22[) only when arg a is close to | . Note that 



1 1 2 • 1 2 2/ 2 2\ 

1 + a sinh r = a (z — z ), 
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so that the assumption |1 + a 2 sinh 2 r\ < c translates to 

|._.i = /0(H~ 1 ) forH<l 
1 01 \0{\a\' 2 ) for \a\ > 1 

To obtain estimates near the turning point, we introduce the functions 

i 

/ 



p{z) : 

By rewriting p(z) in the form 



z - zo 



p{z) 



ay/ Z + Zq 



(z - z a y. 



z 2 -l ' 
we can easily obtain estimates, 

(A 26) \v (k) (z)\ - l^ +k for l a l ^ 1 and l z - z o| =0(|a|- 1 ), 

1 ' ' \\a\ 3+2k for |a| > 1 and \z - zo\ ^ 0(\a\~ 2 ). 

Using the definition of <j> as f s/J dz, we can write q(z) in the form 



t 2 p(z + t(z - z )) dt. 



Then from (j A.26|) we can derive estimates of the same form for q(z) 
(A.27) W k) (z)\ 



|a|t +fc for \a\ < 1 and \z - z Q \ = 0(\a\- r ), 
|a| 3+2fc for \a\ > 1 and \z - z \ = 0(\a\- 2 ). 



Using (|A.26|) and (|A.27|) , with the fact that f /( = _p 2 (|q) 3 and the formula for given 



in (|A.15I) . we obtain 

For \a\ < 1, the result is 
(A.28) 



< 



'OQafi) for \a\ < 1 and \z - z \ = OQa]- 1 ), 
0(\a\ 2 ) for \a\ > 1 and \z - z Q \ = 0{\a\- 2 ). 



'\z-z \<G\a\- 1 

For \a\ > 1 the corresponding estimate is 



(A.29) 



\z-z \<C\a\- 



dz = 0(\a\s), 



dz = O(l). 



Now we can combine the estimates of contributions to ^(r) from all three zones, namely 
(|A~24l) . (|A~25l) . (TOg]) . and (jA~29|) . to obtain 



(A.30) 



|*i(r)|<C(l + |a|-*), 



for arga € [0, ■§] and r G [0, oo), with C independent of both r and a. Applying the 
resulting estimate of \&(r) in (|A.21[) then gives 



3 e 



fe Rc 



This completes the error analysis in the P case. 
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We turn now to the Legendre Q- function and the proof of (|A.12[) . We want the solution 
to be recessive at r — oo, so we set a = in the ansatz (|A.16[) and impose the condition 

(A. 31) ho = 0(r~ 2 ), as r — > oo. 

Using (|A.9|) and (|A.19|I we have 

^ r Ai(fcfO ~ e -(**-*)r 



[1 + a 2 sinh 2 r] 1 " " 2-7T2 
as r — > oo. From the asymptotic (|6.4p we find 

i 

^-i+k^ costir > T (ka + 1) 



1 • ' - 2 '** (S) *(afalir)- 1 (C//)iWb, 



Hence, for the Q-Legendre function we have 

Q ^+- (coshr) r(te + i) 

which proves (|A.12|) . 

To control ho we use the boundary condition (|A.31|) to transform the differential equation 
(|A.17[) for ho into an integral equation, 



h (k,a,r) 
where 



2ire e 



9 fir 1 *) 2 sinh r 1 
K (r, r')V(r') [h (k, a, r') + M(kH(r'))} A ' f" dr', 



#o(r, r') := Ai(fcic(r'))Ai(fc* e - 2 ^CW) - Ai(fcie"^C(^'))Ai(fc^CW)- 
The consequence is that the analog of (|A.21|) for h^ is 



(A.32) 
with 



fee C 4 ^0 



*o(r) 



< Ce 



- fe Rc ■ 



(r) 



sinh r dr' 



Since \&o(»") = Wi(oo) — &i(r), we can simply apply the estimate (|A.30|) from the P case to 
complete the proof. □ 

The first application we need from Proposition IA. ll is a set of good upper bounds. 

Corollary A. 2. Assuming that \ka\ > 1, Re a > 0, and r £ [7*0,7*1], we have the following 
estimates: 



(A.33) 

and 

(A.34) 



Prf +fea (coshr) 



< 



r(fe + i) 



fe Rc[0(a,r)-p(a)] 



Q"i +fca (coshr) 



< Cfcs | a | 2 p _ feRe [^( a>r .)_ g ( a )] 



|r(fca + l)| 



where C depends only on rp and ri . 
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Proof. By conjugation, it suffices to assume that arga € [0, Using the asymptotics, 
(IA.19[) and (|A.20[) . and the first error estimate from (|A.13[) . we have 



(A.35) 



^(^[Ai(^eT() +hi(k,a,r)] 



< Ce 



k Re ( 



> c, for some c > 0, then the estimate (|A.33|) follows 



If we assume that |1 + a 2 sinh 2 r\ 
immediately from (|A.11[) . 

On the other hand, if |1 + a 2 sinh 2 r| < c, then by the assumption that r e [ro,ri], we 
deduce that |a|, |(?!>|, and the ratio + a 2 sinh r] are all O(l). For \cj)\ > fe -1 , we can use 
(|A.35j) to complete the estimate. If \cf>\ < fc _1 , then |Ai(fc 3 e~ C,) + h\{k, a,r)\ is bounded 
by (|A.13[) and the fact that Ai(u>) is regular at the origin. (It is only because of this last 
case that the factor ke must be included in the final estimate.) 

The argument for (|A.34[) is essentially identical. □ 

Our second application of Proposition lA.il is to control the ratio of Legendre functions. 

Corollary A. 3. Assuming that \ka\ > 1, e > 0, and arga 6 [0, 5 — s], we have uniform 
bounds for k sufficiently large: 



Q (coshr) 



Q" i +fcQ (coshr) 



r(fca + l) 



a?r(fc+ 1) 



Ji Rc[2t/)(a,r)— p(a)— q(a)] 



meaning the the ratio of the two sides is bounded above and below by constants depending 
only on e. (The upper bound extends to arga = ^, but the lower bound does not.) 

Proof. By (|A.11|) and (|A.12jl we have 



P 



-k 



(A.36) 



Q' 



_ fcQ (coshr) a hy(k a + l) 
(coshr) _C T(fc + 1) E 



k[p(a)+q(a)] 



Ai(fcle 2 3 S ^) + hi(k,a,r) 
Ai(fciC) +h 2 (k, a.r) 



For some c > 0, consider first the case where \ k<p\ > c. The assumption that arga is bounded 
away from ^ implies that arg£ G [— | — 6%], so that we can apply (|A.19|) to estimate 
both of the Airy functions in (|A.36[) . By choosing c sufficiently large, we can assume that 



the factor l + 0(|u>| 3 / 2 ) appearing in (|A.19[) is bounded away from zero, since 



IM 



in our case. By the estimates (|A.13I) and the assumption \ka\ > 1, by choosing k sufficiently 
large we can assume that \h\\ and \h2\ are arbitrarily small relative to the Airy function 
estimates. Under these assumptions we have 



Ai(fc3 e 2 3 ! () -f hi(k,a,r) 



Ai(fc3() + h 2 (k,a,r) 



o 2fcRc0(o;,r) 



The bound then follows immediately. 

If \kcf>\ < c, then because Ai(iu) is non-zero near the origin, the bound follows immediately 
from (jA.llj) . (|A.12p . and (|A.13|) . provided that k is sufficiently large. □ 
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